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Abstract 

We study a two-dimensional spinless particle in a disordered gaussian 
magnetic field with short time fluctuations, by means of the evolution 
equation for the density matrix < ■x} 1 '\p(t)\x.^ >; in this description 
the two coordinates are associated with the retarded and advanced paths 
respectively. The static part of the vector potential correlator is assumed 
to grow with distance with a power h; the case h — corresponds to a 
(^-correlated magnetic field, and h = 2 to free massless field. The value 
h = 2 separates two different regimes, diffusion and logarithmic growth 
respectively. When h < 2 the baricentric coordinate r = (l/2)(x' 1 '+x' 2 - ) ) 
diffuses with a coefficient D r proportional to x , where x is the relative 
coordinate: x = yS 1 ^ — x*- 2 '. As h > 2 the correlator of the magnetic 
field is a power of distance with positive exponent; then the coefficient D r 

— 2 2 

scales as x . The density matrix is a function of r and ^-,and its width 
in r grows for large times proportionally to log(-^j). 
PACS numbers:71.55.j, 73.20.f 

1 Introduction 

In recent years various studies have been devoted to propagation in the presence 
of a stochastic drift term. The most known example comes from fluid dynamics: 
the transport of a passive scalar in a velocity field. It has been shown that this 
problem can be solved in various physically relevant situations. In particular 
Gawedzky and Kupiainen Q gave a solution in the three dimensional case by 
means of an explicit resummation of the perturbative series. Correlators of 
the scalar field of any order can be obtained. The basic ingredient is the as- 
sumption of zero correlation time for the velocity, which allows for an effective 
evolution generator. Quotation of previous work on the passive scalar can be 
found, e.g., in Rcf. A second, related example is quantum motion in a dis- 
ordered magnetic field. This topic has been raised in various physical contexts: 



Mott-Hubbard systems M, vortex lines in superconductors ||, Quantum 
Hall Effect at filling factor 1/2 In the last case the magnetic field is 

static with zero mean and numerically there is evidence of a derealization tran- 
sition Q,||],[[lO), at odds with a supersymmetric treatment leading to a sigma 
model with unitary symmetry IjOL Fluctuations in the topological density have 
also been called responsible for the derealization transition [|l2| . Recent work, 
in examining the connection between passive scalar and quantum particle in a 
magnetic field, related dclocalizcd states with random antisymmetric disorder, 
but also with symmetric disorder, provided a peculiar sublattice decomposition 
is allowed Quantum mechanics with an imaginary magnetic field and a 

disordered potential, introduced by Hatano and Nelson recently attracted 
much attention, (see, e.g., Efetov fl5fl , showing how nonhcrmiticity can sustain 
extended states in spite of disorder. In a previous work on static magnetic disor- 
der we have shown that it translates, at the semiclassical level, into an imaginary 
(nonlocal in time) magnetic field. We further interpreted the effective action in 
terms of Coulomb gas, and discussed how its strong coupling regime can enhance 
phase coherence among trajectories jL6). Static disorder cannot be studied in 
terms of an effective time evolution generator, due to the non-locality of the 
effective action. Time-dependent disorder, in the limit of zero-correlation time, 
can instead be studied: one can think of a low energy particle moving in a sea 
of highly energetic ones, each carrying a magnetic flux line, or to a particle, 
as, e.g., a vortex, in a fastly varying effective magnetic field. Cooperon dynam- 
ics with time-dependent magnetic fluctuations has been studied by Aronov and 
Wolfle O], in the context of high-T c superconductors. In the present work we 
analizc rapidly fluctuating magnetic disorder on an otherwise free particle: our 
aim is to determine whether the ballistic behavior is frustrated in such a case, 
and what kind of presumably slower propagation sets in. By taking the corre- 
lation length of the fluctuations as the large scale of the system we establish 
the effective (annealed) dynamics for various power-law disorder correlators. In 
the second Section we determine the time evolution generator for the density 
matrix from the Feynman path integral. This is a two-particle description, since 
both retarded and advanced paths must be taken into account. Averaging over 
disorder leads to self-interaction as well as to mutual interaction between the 
paths: the former amounts to a renormalization of the single particle dynamics 
(i.e. it influences the quantum amplitude), while the latter, referring to interfer- 
ence between counterpropagating paths, is essential to quantum probability. In 
the third Section we discuss the particle-antiparticle relative dynamics, which is 
related with the quantum fluctuations in the momentum space. In Section 4 we 
analize the full problem and give the main results on the particle diffusion and 
subdiffusion. In Section 5 we make some comparison with similar behaviors in 
different physical contexts. 



2 



2 Density matrix evolution 



We assume that the correlation time of magnetic fluctuations can be neglected 
and take the following gaussian correlator for the vector potential, in the trans- 
verse gauge: 

D tti/J (k)-5(i) ^< A a Af3 > (M) = ASlpik) — - ^y +h/2 S{t), (1) 

where k^ 1 is the correlation length, h is positive and 5^ a = — (k a kp)/k 2 . 
We will take into account the range < h < 4; the limit h — > , fco — » 
gives a magnetic field B delta-correlated in space, h = 3 can be associated 
with the anomalous skin effect p?|] ; h — 2 gives a free massless -B and, 
apart from logarithmic corrections, corresponds to the enstrophy cascade in 
two-dimensional turbulence (with the vector potential playing the role of ve- 
locity). Notice that as ft. > 2 the correlator of B is a power of distance with 
positive exponent. We write the Feynman path integral representation for the 
product of amplitudes: 

F(x,y;x',y';t) ^< x | t/«(t) |y><x'| U^{t) \ y' >*, (2) 

where EA*'(i),i = 1,2 are the evolution operators of two identical copies of 
single particle systems in the presence of the magnetic field. The time-evolved 
density matrix is given by: 

< x|p(t)|x' > = J dydy'F(x,y;x',y';t) < y|p(0)|y' > (3) 

The representation of Eq. involves a two-particle action: advanced and re- 
tarded path respectively. Averaging over the vector potential couples the paths 
and generates the following effective lagrangian £ e //(X , X) , X = (x^ , x^ 2 )) : 

C ef f = -Xa 3 X + -JL-Xa 3 Da 3 X = -XMX. (4) 

Here g = - , the Pauli matrix 173 acts on the particle indices and the (4X4) matrix 
D = Z)w) (z,j = 1,2) corresponds to the static part of the correlator defined 
in Eq. 0: D^V = = D a>f3 (x.W = ^);D^ = = LWxW - 

x^ 2 )); (a, (3 = x, y). The evolution of the density matrix is given by Eq. [|, where 
F is converted into its average over disorder; in the effective action the imaginary 
part of M is obviously related with dissipation. Since, contrary to the static 
case (see, e.g.,Ref. |lB]), we have locality in time, it is possible to extract a time 
evolution generator H from Eq. |^. This can be performed with no ambiguity 
in a flat metrics, while in general operator ordering prescriptions are needed. 
As discussed at length in Ref . [[l7] , if one evaluates the metrics at the midpoint 
between the initial and final configurations, and consistently normalizes the 
intermediate gaussian integral, one obtains a symmetrized generator. We have: 
H = |(VF- x V) a where F = ±M and (ViV) s = ±(VVi + 2VAV + iW). 
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Here the sum over indices is understood and the gradient is with respect to X. 
The kernel F~ 1 is given by: 

F- 1 = C-^ + i-fm-G- 1 (5) 

m m 

G a ,p = ^ + (^) 2 '[^(0)-<,(xW-xP))] 

This simple form results from D a ^ being diagonal when its argument is equal to 
0. The operator G is diagonal in the particle indices and, as one easily verifies, 
is even under particle exchange; it obviously commutes with D = D(x^ — 
x^ 2 ) ) . By taking into account the transversality condition one ends up with the 
following evolution equation for the density matrix p(x^ , x' 2 ' ; t) in coordinate 

Hp- 

+l(^) 2 Dl, 7 (G-%AdPd? + dfdW) 
+ ^[d a (G- 1 ) a M(df ) +df ) ) (6) 

+ ±(JL)2 d ,[(D -D*) a MG- 1 )^] 

Here the partial derivatives, unless labeled with the particle index, are taken 
with respect to the relative coordinate x = x^ 1 ) — x^ 2 ) and D° , D x refer to the 
static part of the correlator at separation and x respectively. Notice that the 
particle mass is renormalized, and an effective magnetic field (proportional to 
(— ) 2 ) acts with opposite signs on the two trajectories. The remaining terms are 
dissipative, of diffusion, potential and drift type: obviously this acts as an imag- 
inary magnetic field (proportional to — ) acting identically on both trajectories. 
If, in the original hamiltonian, one disregards the A 2 term, one can average the 
two-particle evolution operator emerging from Eq. |^. Following the lines of Ref. 
[Q, one starts with the Neumann series for this operator, taking as perturbation 
the A ■ p coupling, then uses Wick's theorem and averages each term. The 
resulting series can be reexponentiated to an effective two-particle generator, 
which coincides with Eq. |^, but with G = 1. This case is fortunate since the 
transversality condition guarantees that the generator H does not depend on 
ordering prescriptions. 



representation: 
dp 

at 



H = 
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3 Particle-antiparticle relative dynamics 



Upon isolating the dynamics of the relative coordinate, as from the center of 
mass of the system, we study the interaction between the retarded and advanced 
paths. If only the relative coordinate survives, the operator H reduces to: 

+ (9-)*(D° - D*) ^[dJG- 1 ) ^ p ]dp (7) 

+ I ( m )29 ^ (jD ° ~ DX hMG-%^] 

We now show that Eq. |^ in the case of long range fluctuations (feo << 1) 
further simplifies, so that the drift and curvature contributions can be neglected. 
In order to proceed, it is first necessary to examine the disorder correlator in 
that regime. In space coordinates the static part has the form: 

A,,/»(x) = {^){X (x)'8 ai p+X 2 (x)''E a ,f,) (8) 
Ei,i = -Ea,a = cos(2V>) ; Ei, a = E 2 ,i = sin(2^) 

where x and ip are the polar coordinates of x and the JjS are Bessel functions. 
The small ko behavior of Eq. || is readily obtained: 

D °,f> = (^X^W + A^) (9) 

A = d) h R ,Q< h <2 i = k -x 
A = (gff ,2 < h < 4 

87rr(l + /i/2) L V 1 ' J 

The explicit form of G is then; 

G = ^-{ J ^^?A^ + A) (10) 



The subtracted correlator has a power-law behavior, and is 0(fc§) in the limit 
fco — > when h < 2: 

D°-D X = -^A (11) 
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The ratio j± has the dimensions of a momentum, let us call it II. The order of 

fc 

magnitude of the A 2 contribution can be estimated by means of the adimen- 
sional coupling constant 2 - J jp-, where r e is the classical radius of the electron 

(r e = We define the large fluctuations regime, where in G we can disregard 
the identity with respect to the second term: 

g 2 ■ A _ r e ■ n 



■ h ■ k, 



» 1 (12) 



Notice that this relates the amplitude of the fluctuations (associated with A) 
with the correlation length fc^ 1 . The kernelF -1 , a function of II, h and of the 
quantum flux unit ^ = (f^y, reduces to: 

1 
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2-n-h 



It is easily verified that in both cases (h < 2 , h > 2) the matrix A satisfies 
the transversality condition; this implies that the generator is, independently of 
ordering prescriptions: 

6 = ~ \ A ^) d « ^ ( 14 ) 

2 g 11 • 7 7 

The dominant contribution is pure diffusion with a coefficient 0(fcg ), while the 
first correction has scaling form, and in both cases it corresponds to a faster 
spreading. When < h < 2 , the correction is of order k 2h and, taken alone, 
would give superdiffusion (t w x 2 ~' 1 ); when 2 < h < 4, it has order kQ +h 
with t w (In a;) 2 . A complete analysis of Eq. |ljis needed in order to understand 
the interplay between diffusion and this faster mechanism; we leave it to future 
work. We finally discuss the regime of weak fluctuations ( 2 - J jp- << 1)- When 
in G we disregard the second term, F^ 1 has a scaling form: 

p- 1 = -(£) 2 ni (15) 

Recall that while in the h < 2 range this is 0(/cq), in the complementary range 
it diverges as k^~ h in the limit fco — > 0. Corrections in the small fluctuation 
parameter come from the expansion of G _1 . Since only the dominant term 
preserves the transversality condition, the operator H has a drift as well as a 
curvature term. The dominant part gives rise to superdiffusion (x w t 1 /( 2_/l )) 
when < h < 2 and time exponential behavior (as previously illustrated) 
when 2 < h < 4. This amounts to disregarding the A 2 term . The generator 
corresponds then to the one obtained for the second order correlator of a passive 
scalar by Gawedzki and Kupiainen JjJ . We stress in conclusion that the scaling 
form of the disorder does not transfer to the generator, due to the fluctuations of 
A 2 ; this reduces the strong tendency of the two trajectories towards separation, 
from fluctuation-supported superdiffusion to diffusion. 
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4 General case 

Let us write equation |^ in terms of the relative coordinate and of the baricentric 
coordinate r = ^(x' 1 ) +x^ 2 ^) : 

I - *' <«> 

H = 7^—( G ~)a,l3(-^ p h T. q ) + [fl ( G ~ )<*>P\a 

2m Or a ox p Ox a orp 2m ox a orp 

Here the limit fco — * gives a finite nonzero result provided that h < 2: as 
discussed in the previous section, the diffusion term for x is 0(fcg ), similarly the 
other terms go to zero in the limit. From Eq. [tfj we get then: 

If one disregards the anisotropy, one easily concludes that this implies an evo- 
lution through the combination r V" . Let us consider: 

H Pa V(x/2)- h V 2 r (18) 

where I? is constant. Let us take as initial condition a gaussian in both particle 
coordinates with equal width: 



p(r, x) = (^-) 2 exp(~w(r 2 + (x/2) 2 )) (19) 



The solution is: 



P(r ' x ' t} = 2p l/2«, + 2P(2/*)^- (w/4)a;2 - (1/2) (l/2«,) + 22?(2/x)"t] 

(20) 

The Wigner function, which in the classical limit goes into the phase space 
distribution, can be evaluated at zero momentum; it is defined as: 

W(r, k,t) = -^J d 2 xexp{-ik ■ k) P (v, x, t) (21) 

In the large time limit we disregard l/(2u>) with respect to 2T>(2/x) h t and 
obtain: 

I poo 

W(r,k = 0,t) = - , 2oh l" / dx-x h+1 exp(~(w/4)x 2 -Qx h ) 



(2tt) 2 2 h + 2 Vt J 



w 1 rW 1 " 



(2tt) 2 2 h + 2 Vt J 



dxx h+1 ex P (-Qx h ) 



= r 2 ■ (2 hJr2 VtY x Vt > - • — < 2 h ■ w h/2 
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The integral can be explicitly written in terms of the truncated gamma function, 

2 

and turns out to be a function of ^r. The k dependence of W will not be 
discussed in detail here. Notice only that, for realp, exp(—a - \p\ h ) is the Fourier 
transform of the Levy distribution Lh(z) ks z(h 1 +1) (z » 1) Q. The density 
matrix displays this stretched exponential for every initial condition provided 
it is gaussian-shaped in r. It is reasonable to assume that the large k behavior 
of W is only marginally influenced by the initial x dependence (x << 1), if 
sufficiently smooth. The Wigner function reduces then to the Q-derivative of the 
Levy distribution in momentum space , where the derivative absorbs the factor 
x , from the normalization in r . Thus the relative motion can be depicted as 
a sequence of jumps and stopovers in momenta, giving rise to a hyerarchy of 
clusters. So far we have been dealing with the fco - > limit; as ko is small but 
different from zero we must take into account other terms, starting with the 
first correction, which gives pure diffusion in x (see the previous section). We 
discuss a situation of this sort in the regime 2 < h < 4. Now D r is 0(k,Q~ 2 ); 

,h 

this, together with the next to leading term gives ( << 1): 

6 ~ { ~A • K" 1 ^) (^(r- 1 )-,^^ + (^K] (22) 

Disregarding anisotropy and Fourier transforming with respect to r we get a 
particle in a centrifugal barrier: 



% = [-(f) 2 + ^]p (23) 

r = (V^. tj 
9 A 

q 2 = (4/3)(/i-2)^, 

where K is the momentum with respect to r. Let us search for a similarity 
solution, p — f(u,q) , u = — , as suggested by the scaling form of the equation; 
the result is: 

f" + (- + 1/4)/' - (f )/ = (24) 

Two independent solutions (when q is not an integer) are: f(u, q) = u s F(s, 2s + 
I; — u/4), where s = q/2 , — q/2 and F is the confluent hypergeometric function. 
In order to avoid singular behavior at small u we exclude the negative s case. 
We numerically computed the Fourier transform , in order to reconstruct the 
r dependence; in Fig.l we show p(r , u) for a set of values u < 1, i.e. in the 
large time regime. We find that the width in r increases, as u approaches zero, 
linearly in the logarithm of u^ 1 , up to four orders of magnitude (see Fig. 2). 
Disregarding the next to leading term would result, similarly to the previous 

2 2 

case, in the combination r ' x . 
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5 Conclusions 



We studied the two-dimensional dynamics of a particle in the presence of a 
stochastic magnetic field in the fast fluctuation regime. We have found a transi- 
tion by varying the exponent of the disorder correlator. Disregarding the bari- 
centric motion, we find that diffusion is the dominant behavior in the ko < < f 
regime. Superdiffusion arises in the first correction, with the scaling t w x 2 ~ h 
when h < 2 and t « (In a;) 2 when 2 < h < 4. The dominance of diffu- 
sion originates in the A 2 term; if this term is neglected we obtain fluctuation- 
sustained superdiffusion with the previously mentioned power laws. In other 
words the tendency of the retarded and advanced trajectories to spread apart 
very rapidly (favored by the A • p term) is partly stabilized by the quadratic 
term. In the general case, the limit fco — > can be performed exactly in the range 
h < 2. The solution (see Eq. [|(]) shows diffusion in r, with the dependence 
r ~ (t/x fl ) 1 ^ 2 . The Wigner function W at zero momentum, i.e. the average 
over particle separation, describes diffusion (t « r 2 ); at k large it appears in 

2 

the form W ~ Wy^r-^K), its k behavior being associated with isotropic Levy 
flights in the momentum space. Random walk in a quantum mechanical system 
was found in the Harper model at its critical point Jl9|,p0||. It was also de- 
rived from a two-dimensional anisotropic random lattice model, describing edge 
propagation along the surface of a quantum Hall multilayer |2l| , |2^| , |p3| , Q . 
In that case backscattering is allowed in the direction of the magnetic field and 
neglected in the orthogonal direction; when the latter is interpreted as time one 
has a one-dimensional chain with time-dependent hopping disorder where the 
quantum particle undergoes diffusion. As long as hopping mimicks the presence 
of a magnetic field, this model can be taken as a one-dimensional version of 
ours. Propagation is further inhibited when 2 < ft.; by including the leading 
and next to leading terms we determined the Fourier transform with respect 
to r of the density matrix (see Eq. and arguments following it). From the 
large time behavior of p(r,x, t) one extracts r w log(-p): this is a much slower 
spread than expected from the leading term alone, giving r « (t/x 2 ) 1 ^ 2 . This 
ultraslow, logarithmic diffusion occurs in biased random motion on percolation 
clusters or globally isotropic fractals (see and references therein) . Quantum 
mechanically, it was found by suitably perturbing at a single point an Ander- 
son hamiltonian in the localized phase Q. When disorder is static, p7|,p8[, 
smoothly varying intense magnetic fields confine nearly free states within a nar- 
row one-dimensional region along the zero field lines, such lines thus providing 
a quantum percolation network. The present case can be more properly de- 
picted in terms of abrupt jumps in the space of momenta, with loss of phase 
coherence. This apparently results in random walk when the correlator of the 
magnetic field B decays with distance {h < 2) and in logarithmic behavior 
when it grows (2 < h). 
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Figure captions 

Fig.l: Plot of the solution p(r,u) (see Eqj2^) at different values of u as 

approaches zero (large time behavior). 

Fig. 2: Width in r of the solution versus log(u). 



References 

fll K.Gawedzki,A,Kupiainen: Universality in turbulence: an exactly soluble 
model. Lecture Notes of the 34 th School of Theoretical Physics in Schlad- 
ming,Austria, March 4-11,1995. 

G.Falkovich,I.Kolokolov,V.Lebedev,A.A.Migdal,Phys.Rev.E54, 
4896(1996). 

L.B.Ioffe,A.I.Larkin, Phys.Rev.B 39,8988(1989). 
L.B.Ioffe,.V.Kalmeyer, Phys.Rev.B 44,750(1991). 
N.Nagaosa,P.A.Lee,Phys.Rev.Lett. 64,2450(1990). 

M.V.Feigelman,Physica A 168, 319(1990). 

E.M.Chudnovsky, Phys.Rev. 51,15351 (1995). 

V.Kalmeyer, S.C.Zhang, Phys.Rev.B 46,9889(1992). 

B.I.Halperin,P.A.Lee,N.Reed,Phys.Rev. b 47,7312(1993). 

Y.Avishai,Y.Hatsugai,M.Kohmoto,Phys.Rev. b47,9561(1993). 

V.Kalmeyer,D.Wei,D.P.Arovas,S.C.Zhang,Phys.Rev. B 48, 11095(1993). 

D.Z.Liu,X.C.Xie,S.Das Sarma,S. C.Zhang, Phys.Rev. B52, 5858(1995). 

A.G.Aronov,A.D.Mirlin,P. Woffle, Phys.Rev.B 49,16609(1994). 

S.C.Zhang,D.Arovas, Phys.Rev.Lett. 72, 1886(1994). 

J.Miller,J.Wang, Phys.Rev.Lett. 76, 1461(1996). 

N.Hatano,D.Nelson, Phys.Rev.Lett. 77,570(1996). 



K.B.Efctov, |c;ond-mat/9706055| , June 6th,1997. 
B.Cardinetti,V.G.Benza,L.Molinari, Phys.Rev. B52, R5499(1995). 
A.G.Aronov,P.Wolfle,Phys.Rev. B 50,16574(1994). 



10 



[18] E.W.Montroll.B.J.West, in:" Fluctuation Phenomena", 

E.W.Montroll,J.L.Lcbowitz eds., North Holland, Amsterdam, 1987. 

[19] H.Hiramoto,S.Abe,J.Phys.Soc.Jpn. 57.230(1988). 

[20] T.Gciscl,R.Kctzmcrick,G.Petschel, Phys.Rcv.Lett. 66,1651(1991). 

[21] J.T.Chalker,A.Dohmen,Phys.Rev.Lett. 75,4496(1995). 

[22] L.Balents,M.P.A.Fishcr,Phys.Rev.Lett. 76,2782(1996). 

[23] L.Balents,M.P.A.Fisher,M.R.Zirnbaucr, Nucl.Phys. B483,601(1997). 

[24] H.Mathur,Phys.Rcv.Lctt. 78, 2429(1997). 

[25] J.P.Bouchaud,A.Georges,Phys.Rep. 195,127(1990). 

[26] R.del Rio, S.Jitomirskaya,Y.Last,B. Simon, Phys.Rev.Lett. 75,117(1995). 

[27] J.E.Muller, Phys.Rev.Lett. 68,385(1992). 

[28] D.B.Chklovskii,P.A.Lec, Phys.Rev. B 46,18060(1993). 



11 




12 



